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The optimal redshift for detecting ionized bubbles in HI 21-cm maps 
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ABSTRACT 

The detection of individual ionized bubbles in HI 21-cm maps is one of the most promising, 
direct probes of the epoch of reionization (EoR). At least 1000 hrs of observation would 
be required for such a detection with either the currently functioning GMRT or the upcoming 
MWA. Considering the large investment of telescope time it is essential to identify the "optimal 
redshift" where the prospects of a detection are most favourable. We find that the optimal 
redshift is determined by a combination of instrument dependent factors and the evolution of 
the neutral fraction a;Hi- We find that the redshift range 8.1 ± 1.1 and 9.8 ± 1 are optimum 
for detecting ionized bubbles with the GMRT and MWA respectively. The prospects of a 
detection, we find, are more favourable in a scenario with late reionization with xhi ~ 0.5 at 
z ~ 7.5 as compared to an early reionization model where xhi ~ 0.5 at z « 10. In the late 
reionization scenario, for both instruments a Scr detection is possible for bubbles of comoving 
radius Rb > 30 Mpc with 1000 hrs of observation. Future observations will either lead to 
the detection of ionized bubbles, or in the event of non-detection, lead to constraints on the 
product .thi R2 for the observational volume, where 7 = 1.5 and 2 for GMRT and MWA 
respectively. 

Key words: cosmology: theory, cosmology: diffuse radiation. Methods: data analysis 



1 INTRODUCTION 

It is currently accepted that the Universe was reionized by the 
growth of ionized bubbles around luminous sources in the redshift 
range 2; ~ 6 — 15 (Fan et al. 2006; Choudhury & Ferrara 2006; 
Komatsu et al. 2008). Detection of individual ionized bubbles (HII 
regions) in HI 21-cm maps of reionization is one of the major, im- 
portant approaches that will be adopted by the present and upcom- 
ing radio experiments (GMRT, MWA, LOFAR, SKA) to probe the 
EoR. Such observations will directly probe the properties of the 
ionizing sources and the evolution of the surrounding IGM (Wyithe 
& Loeb 2004; Wyithe, Loeb & Barnes 2005; Maselli et al. 2007; 
Geil & Wyithe 2007) and are expected to complement the study of 
reionization through the power spectrum of HI brightness tempera- 
ture fluctuations. Detection of individual bubbles is a big challenge 
because the HI signal will be buried in strong foregrounds and sys- 
tem noise (Ali et al. 2008). 

In an earlier paper (Datta, Bharadwaj & Choudhury 2007), 
hereafter referred to as Paper I, we have proposed a visibility based 
matched filter technique to optimally combine the entire HI sig- 
nal from an ionized bubble while removing the foregrounds and 
minimizing system noise. Using visibilities has an advantage over 
image based techniques because the system noise contribution in 
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different visibilities is independent whereas the noise in different 
pixels of a radio-interferometric images is not. Our investigations 
show that for both the GMRT and the MWA, at redshift z = 8.5, 
it will be possible to detect ionized bubbles of comoving radius 
Kb > 40 Mpc and Rt > 22 Mpc in 100 and 1, 000 hours of obser- 
vations respectively. We also find that fluctuations in the HI outside 
the bubble that we are trying to detect impose a fundamental restric- 
tions on the smallest bubble that can be detected. Assuming that the 
HI outside the bubble traces the dark matter, we find that it will not 
be possible to detect bubbles with comoving radius less than 8 and 
16 Mpc with the GMRT and the MWA respectively, however large 
be the integration time. In a subsequent paper (Datta et al. 2008), 
hereafter referred to as Paper II) we have used simulations to val- 
idate our matched filter technique and assess the impact of patchy 
reionization outside the bubble that we are trying to detect on the 
bubble detection. 

The question "What is the optimal redshift for bubble detec- 
tion?" is particularly important when planning future observations. 
Estimates show (Papers I and II) that at least 1000 hrs of obser- 
vation will be required for a detection with either the GMRT or 
the MWA. Considering the large investment in observing time, it is 
important to target the redshift where the prospect of a detection is 
most favourable. In addition, it is important to have a clear picture 
of the different factors that contribute towards deciding the most 
optimal redshift. We expect this to provide insights useful for the 
design of future observational programmes and also the design of 
future low-frequency radio telescopes. 
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A variety of redshift dependent factors influence tiie signal 
from an ionized bubble. While a number of these pertain to the 
instrument in question, the frequency dependence of the sky tem- 
perature and the redshift evolution of the neutral hydrogen fraction 
also play an important role. In this paper we analyze all of the ef- 
fects that determine the optimal redshift for detecting ionized bub- 
bles. We consider two different models for the redshift evolution 
of the neutral fraction and make predictions for the GMRT and the 
MWA. 

The paper is organized as follows. In Section 2 wc briefly re- 
view the matched filter technique for bubble detection. In Section 3 
we establish scaling relations for the matched filter signal to noise 
ratio (SNR) assuming an uniform baseline distribution. We also dis- 
cuss the models of HI evolution that we adopt. We present our re- 
sults and conclusions in Section 4. 

Throughout out this paper we adopt cosmological parameters 
from (Dunkley et al. 2009). For the GMRT we use the antenna spec- 
ifications from their web site and for the MWA we use the instru- 
mental parameters from (Bowman et al. 2007). 

2 THE MATCHED FU^TER TECHNIQUE FOR 

DETECTING IONIZED BUBBLES IN REDSHIFTED 
21-CM MAPS 

The visibility recorded in a radio-interferometric observation of an 
ionized bubble can be written as 

V{U, v) = S{U, v) + HFiU, v) + N{U, v) + F{{j , v) . (1) 

Here we refer to C/ = d/A as a baseline, d being the physical 
separation between a pair of antennas projected on the plane per- 
pendicular to the line of sight and A is wavelength corresponding 
to the observed frequency v. In eq. (1) Sijj ,v) is the HI signal 
from the ionized bubble, HF{U, v) is the contribution from fluc- 
tuations in the HI outside the target bubble, N{U, u) is the sys- 
tem noise and F{U , v) is the contribution from other astrophysical 
foregrounds. The contributions HF{U, v),N{U, v) and F{U ^ v) 
are all assumed to be random variables with zero mean, whereby 
{V{U = {S(U The angular brackets here denote average 
with respect to different realizations of the HI fluctuations, system 
noise and foregrounds. 

We consider a spherical ionized bubble of comoving radius 
Rb centered at redshift Zc located at the center of the field of view 
(FoV). The bubble is assumed to be embedded in an uniform IGM 
with neutral hydrogen fraction xhi- A bubble of comoving radius 
Rb will be seen as a circular disc in each of the frequency chan- 
nels that cut through the bubble. At a frequency channel v, the 
angular radius of the disc is 6,, = {Rb/ri,)\/l — {Au/Aub)^ 
where Av = fc — is the distance from the bubble's center 
= 1420MHz/(l -I- Zc) and Aub = Rb/rl is the bubble's 
radius in frequency space. Here r,, is the comoving distance cor- 
responding to 2: = (1420MHz/i^) — 1, and = dvu/dv. The 
expected visibility signal S{U , v) in each frequency channel is the 
Fourier transform of a circular disc which can be expressed in terms 
of Ji{2tvU6v) , the first order Bessel function (Paper I). In each 
channel, the signal has a peak value | 5(0, i^) |= irxniTiyOt where 
/. = 2.5 X 10^^ (M) (_||.) ,s the background HI 

specific intensity expected from completely neutral medium. The 
signal is largely contained within baselines U < Uo = 0.61/01/ 
where the Bessel function has its first zero crossing, and the signal 
is much smaller at larger baselines. The signal S{U,v) picks up an 
extra phase if the bubble is shifted from the center of the FoV. The 
amplitude of the signal also falls because of the telescope's primary 



beam pattern (Paper I), and in this paper we restrict our analysis to 
the most favourable situation where the bubble is at the center of 
the FoV. The terms xm, Tv, 0^ and Aub are all redshift dependent, 
and hence the signal too is strongly redshift dependent. In Section 
3 we will discuss the combined effect of all these factors on bubble 
detection. 

In order to detect an ionized bubble whose expected signal is 
S{U, v) we use the matched filter Sf{U, v) defined as 



e i-2i 



'vc-B'/2 



(2) 



Note that the filter is constructed using the signal that we are try- 
ing to detect. The term (u/iyc)'^ accounts the frequency dependent 
U distribution for a given array. The function G is the Heaviside 
step function. The second term in the square brackets serves to re- 
move the foregrounds within the frequency range i/c — B'/2 to 
i^c + B' /2. Here B' = 4 Avb is the frequency width that we use to 
estimate and subtract out a frequency independent foreground con- 
tribution. This, we have seen in Paper I, is adequate to remove the 
foregrounds such that the residuals are considerably smaller than 
the signal. Further we have assumed that B' is smaller than the 
total observational bandwidth B. The filter Sf{U, u) depends on 
[Rb, Zc, dc] the comoving radius, redshift and angular position of 
the target bubble that we are trying to detect. 

Bubble detection is carried out by combining the entire ob- 
served visibility signal weighed with the filter. The estimator E is 
defined as 



E = 



J2S}0a,l^b)V{Ua,l^b) 



El 

a,b 



(3) 



where the sum is over all frequency channels and baselines. The ex- 
pectation value (E) is non-zero only if an ionized bubble is present, 
and it is zero if there is no bubble in the FoV. 

The system noise (NS), HI fluctuations (HF) and the fore- 
grounds (FG) all contribute to the variance of the estimator 



({AEr) = [iAEf)^^ + (^iAE) 



(AEf) 



FG 



(4) 



A 3cr detection is possible only if (E) > 3y ((Ai?)^). In a situ- 
ation where this condition is satisfied, the observed value Eq may 



be interpreted as a detection if Eq > Sy {{AE)^). 

Because of our choice of the matched filter, the contribution 
from the residuals after foreground subtraction ( (AE)'^ ) is pre- 

\ / FG 

dieted to be smaller than the signal (Paper I) and we do not consider 
it in the subsequent analysis. The contribution ( (AE)^ ) which 

\ / HF 

arises from the HI fluctuations outside the target bubble imposes a 
fundamental restriction on bubble detection. It is not possible to de- 
tect an ionized bubbles for which (E) < ^ ^(AiJ)^^ . Bubble 

detection is meaningful only in situations where the contribution 

from HI fluctuations is considerably smaller than the expected sig- 
nal. Once this condition is satisfied, it is the SNR defined as 



SNR = {E}/^J{{AE)'^}t,s (5) 
which is important for bubble detection. The value of SNR peaks 
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when the parameters of the fiher exactly match the bubble that is ac- 
tually present in the observation, and decreases from its peak value 
if there is a mis-match (Paper II). In the subsequent analysis we 
shall use this to assess the redshift that is optimal for bubble detec- 
tion. 

It is possible to analytically estimate {E), {{AE)'^)t^s and 
((A_B)'^)hf in the continuum limit (Paper I). We have 



and 

[AEf 



dupN{U,u) Sf*{U,v)S{U,v) 



(6) 



dupM{U,u) I Sf{U,u) |2 . (7) 



: j cFU j dv\ j du2 



j dB^^ \ I dB^.^ 



UF J J J \ dT J \ dT 

X Pn{U, Vi)Pn{U, V2)Sf''{ij, Vi_)Sf{U, V2 



(8) 



where By is the specific intensity of blackbody radiation (which 
can be approximated as IksTv^ ((? in the Rayleigh- Jeans regime) 
and Pn{U , is the normalized baseline distribution function de- 
fined so that J d U J dvpN{U ,v) = 1. For a given observa- 
tion, d^U dv Pn{U, v) is the fraction of visibilities in the interval 
d^U dv of baselines and frequency channels. Further, we expect 
Pn{U, v) oc v^'^ for an uniform distribution of the antenna sepa- 
rations d. 

The term o in eq. (7) is the rms. noise expected in an image 
made using the radio-interferometric observation being analyzed. 
Assuming observations at two polarizations, we have 

a = ^"^ ^^""^ (9) 

A^ff^/NbtabsB 

where fes is the Boltzmann constant, Tays the system temperature, 
Aefj the effective collecting area of an individual antenna in the 
array, Ni, the number of baselines, toba the total observing time 
and B the observing bandwidth. 

The contribution from HI fluctuations {(A_B)^)hp is cal- 
culated using (^^), the conversion factor from temperature to 
specific intensity at frequency v, and C2ttu{i'i,V2) the multi- 
frequency angular power spectrum (MAPS; (Datta, Choudhury & 
Bharadwaj 2007)). The HI distribution during the epoch of reion- 
ization is highly uncertain. The value of ((A_E)^)hf is sensitive 
to the size and clustering of the ionized patches outside the target 
bubble (Paper II). Given the lack of information, we make the sim- 
plifying assumption that the HI outside the target bubble exactly 
traces the dark matter. This gives the most optimistic constraints on 
bubble detection, the constraints are more severe if patchy reion- 
ization is included. 

3 SCALING RELATIONS 

The scaling of the expectation value of the estimator with various 
parameters can be estimated from eq. (6) whereby 



{E) oc UlAubu:'' I S{Q,u,) 



(10) 



Here we have assumed that B' is larger than the frequency extent of 
the bubble Ai^b and that the baselines in the array extend well be- 
yond Uo. Further it is assumed that the array configuration is such 
that the antenna separations d are uniformly sampled, whereby 
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Figure 1. The evolution of the mean neutral fraction xhi with redshift for 
the two different reionization models discussed in the text. 



follows from eq.(7) that {{AEf) oc (j^{E). We use these and 
the relations from the previous Section to determine that the SNR 
scales as 



SNRoc AeffVNbtobs— — xm 

J- svs 



Rl (l + Zc 



(11) 



This completely quantifies the dependence on the telescope pa- 
rameters, observation time, system temperature, neutral fraction, 
bubble radius and the background expansion history. In principle, 
measurements of the SNR will provide an unique and independent 
way to probe the source properties ( through Rb\ (Yu 2005) ), inter 
galactic medium, and the background cosmology during the EoR. 
For the redshift range of our interest it is reasonable to assume 
r^, oc (1 + r^, oc (1 + zf'= and H{z) oc (1 + zf"^. Fur- 

ther, for the frequency range of our interest the system temperature 
is dominated by the sky temperature which scales as Tsky oc 



with /3 ~ 2.6 which implies Tg,, 



(1 + z)^ . The effective col- 



lecting area is nearly constant for dish antennas like the GMRT 
whereas it scales as A^ff oc v~'^ for dipoles (eg. MWA). Combin- 
ing all of these factors we determine the scaling of the SNR with 
redshift 



SNRoca;Hi(2)(l + z)" 



(12) 



p{U,v) 



Considering the noise contribution next, it also 



where a = —fi — 1 or a = — /? + 1 for dish antennas or dipoles 
respectively. While xm increases with z, the other term (1 + 2)" 
has the opposite behaviour These two competing effects decide 
the redshift where the SNR peaks which is the optimal redshift for 
bubble detection. 

The baseline distribution, in general, does not uniformly sam- 
ple all baselines. Typically, the sampling falls at larger baselines 
and we do not expect the scaling relations discussed here to be ex- 
actly valid. The deviations from the scaling relations depend on the 
bubble size and the array configuration, and in the next Section we 
discuss these for the GMRT and the MWA. 

3.1 Evolution of neutral fraction with redshift 

In this work, we consider two physically motivated models of 
reionization, namely, the early reionization (ER) and the late reion- 
ization (LR) scenario. These models are constructed using the 
semi-analytical formalism (Choudhury & Ferrara 2005; Choud- 
hury & Ferrara 2006a) which implements most of the relevant 
physics governing the thermal and ionization history of the IGM, 
such as the inhomogeneous IGM density distribution, three differ- 
ent classes of ionizing photon sources (massive Pop III stars. Pop 
II stars and QSOs), radiative feedback inhibiting star formation in 
low-mass galaxies and chemical feedback for transition from Pop 
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Redshift z 

Figure 2. Assuming xjji = 1, the dashed hnes show the predicted scal- 
ing of the SNR for uniform baseline coverage (eq. 12). the solid lines are 
calculated numerically incoiporating non-uniform baseline coverage. For 
both GMRT (upper panel) and MWA (lower panel), the upper curves are 
for Rf, = 50 Mpc with 1000 hrs observation, and the lower curves for 
Rt = 20 Mpc with 4000 hrs. 

Ill to Pop II stars. The models are consistent with various observa- 
tional data, namely, the redshift evolution of Lyman-limit absorp- 
tion systems (Storrie-Lombardi et al. 1994), the Gunn-Peterson ef- 
fect (Songaila 2004), electron scattering optical depths (Kogut et 
al. 2003), temperature of the IGM (Schaye et al. 1999) and cos- 
mic star formation history (Nagamine et al. 2005). We assume that 
these two models "bracket" the range of models which are consis- 
tent with available data. 

In ER scenario, hydrogen reionization starts around 2 « 16 
driven by metal-free (Pop III) stars, and it is 50% complete by z ~ 
10. The contribution of Pop III stars decrease below this redshift 
because of the combined action of radiative and chemical feedback. 
As a result, reionization is extended considerably completing only 
at 2; ~ 6 (Figure 1). In LR scenario, the contribution from the 
metal-free stars is ignored, which makes reionization start much 
later and is only 50% complete only around z ~ 7.5. The main 
difference between the ER and LR models is in their predictions 
for the electron scattering optical depth (which is 0.12 and 0.06 for 
the ER and LR scenarios, respectively). 

4 RESULTS AND CONCLUSIONS 

We consider two possible definitions of the 'optimal redshift' for 
bubble detection. The first is the redshift where, for a fixed observ- 
ing time and bubble radius Rt, the SNR is maximum. Another pos- 
sibility is, for a fixed observing time and SNR, the redshift where 
a bubble of the smallest size can be detected. While the two defini- 
tions are the same if the instrument has uniform baseline coverage, 
we do not expect this to be true in general. 

We have used equations (6), (7) and (8) to calculate the SNR 
and determine the constraints from HI fluctuations. The baseline 
distribution function Pn{U, v), which we assume to be circularly 
symmetric (pn{U, v) = Pn{U, u)), has been calculated in Paper I 
for both the GMRT and the MWA. In both cases pn{U, u) falls off 




6 7 8 9 10 11 12 13 14 15 16 



Redshift z 

Figure 3. The SNR for Rt = 50 Mpc and 1000 hrs observation. Results 
are shown for both GMRT and MWA using the two different reionization 
models (ER and LR) discussed in the text. 
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Figure 4. SNR contours as a function of the redshift z and comoving bub- 
ble radius Ri,, considering 1000 hrs of observation with the GMRT (upper 
panels) and MWA (lower). The left and right panels show the ER and the 
LR scenarios respectively. The shaded region is ruled out due to the HI 
fluctuations. 

with increasing U. For the GMRT Pn(U, v) is roughly constant for 
antenna separations d < 1 km and it extends out to large baselines 
d ~ 25 km. For the MWA we have assumed that the antennas are 
distributed over a circular region of diameter 1.5 km, with the num- 
ber density of antennas falling as 1 /r^ with the distance from the 
center. 

We first consider, for a fixed bubble radius and observing time, 
how the SNR varies with z. Assuming xm — 1 and uniform base- 
line coverage, we expect that SNRoc (1 + z)°' with a = —3.6 
and -1.6 for GMRT and MWA respectively. For the GMRT, we 
find (Figure 2) that the predicted scaling holds for large bubbles 
Rb > 50 Mpc where the entire signal lies within a small baseline 
range which is nearly uniformly sampled. For smaller bubbles a 
significant amount of signal spreads over to larger baselines which 
are not uniformly sampled. We find that a changes, approximately 
linearly, from —3.6 to —4.1 as Rb is varied from 50 to 20 Mpc. 
For the MWA, the non-uniform baseline coverage makes the scal- 
ing steeper than —1.6 for all values of Rb, and we find a — —2.4 
and —2.5 for Rb — 50 and 20 Mpc respectively. We combine these 
findings with earlier results at z = 8.3 (Paper I) as to how the SNR 
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Table 1. Dimensionless parameters required to calculate the SNR using 
eq. (13). The values of Rh are restricted to the range 50 Mpc > Rb > 

20 Mpc. 





K 


a 


7 


GMRT 


9.1 


-3-6 -(50- 5^) /60 


1.5 


MWA 


13.4 


-2.4 


2.0 



Table 2. For R^ = 50 Mpc and 1000 hrs of observation, the optimal 
redshift Zo where the SNR peaks, the peak value and the z range where the 
SNR is within 80 % of the peak value. 







Zo 


Peak SNR 


80 % z range 


GMRT 


ER 


9.2 


3.6 


7-12 




LR 


7.6 


8.4 


6.8 - 9.2 


MWA 


ER 


11.0 


6.59 


8.8- 14 




LR 


8.4 


11 


7.1 - 10.8 



scales with Rb to obtain 

«^^ = -'^ (look)" (^)"(50^)' 

where t is the observing time, and a, 7 and K are parameters whose 
values are Usted in Table 1. This expression is found to match the 
numerically computed SNR to within 20%, which is quite adequate 
given the large uncertainty in xm- 

Considering Figure 3 which shows the SNR for the two reion- 
ization models , we find that it increases monotonically as z de- 
creases when xhi ^ 1 and thereafter declines rapidly once xhi < 
0.5. The peak SNR, the corresponding optimal redshift Zo and the z 
range where the SNR is within 80 % of the peak value are tabulated 
in Table 2. Results have been shown only for Rb = 50 Mpc and 
t = 1000 hrs of observation, these can be easily scaled to other Rb 
and t values using eq. (13). The z dependence is not very different 
for smaller bubbles in the range 50 > J?b > 20 Mpc. 

The effective collecting area of the individual MWA antennas 
increases with wavelength as A^. This reduces the noise at higher 
redshifts, and puts the MWA at an advantage over the GMRT in 
detecting bubbles at high redshifts. This also pushes the optimal 
redshift for MWA to a higher value as compared to the GMRT (Ta- 
ble 2). The MWA is also at an advantage over the GMRT in de- 
tecting large bubbles (Rb ~ 50 Mpc, Figure 3). The SNR scales 
differently with Rb for the two instruments (Table 1), and the ad- 
vantage that the MWA has for large bubbles balances out as the 
bubble size is reduced. GMRT and MWA have nearly comparable 
SNR for Rb = 30 Mpc. 

We next consider the other definition of the optimal redshift 
where for a fixed observing time and SNR, we determine Zo where 
a bubble of the smallest size can be detected. Considering the con- 
stant SNR contours in Figure 4, we find that the Zo values are 
roughly consistent with those in Table 2. This shows that for both 
the GMRT and the MWA, for 50 > Rb > 20 Mpc the two def- 
initions predict the same optimal redshift which is approximately 
independent of the bubble size. We do not expect this to hold for 
smaller bubbles Rb ~ 10 Mpc where a detection is possible only 
with the GMRT, the signal being smaller than the HI fluctuations in 
the MWA (Paper I). 



Given the lack of knowledge about the reionization history, it 
would be most judicious to choose a redshift where a high SNR is 
predicted for both the ER and LR models. We find that the redshift 
range 7 — 9.2 and 8.8 — 10.8 are most appropriate for the GMRT 
and MWA respectively. For both instruments, the prospects of a de- 
tection are considerably improved in the late reionization scenario. 
Assuming 1000 hrs of observation, in the ER and LR models re- 
spectively, a 3cr detection is possible with the GMRT for Rb ~ 50 
and 30 Mpc or larger. The same figures are 40 and 30 Mpc for the 
MWA. 

The actual distribution of bubble sizes is an important is- 
sue for bubble detection. This depends on the reionization his- 
tory and the distribution of ionizing sources which are largely un- 
known. We generally expect a predominance of larger bubbles at 
lower redshifts. Analytic estimates (Furlanetto et al. 2006; Rhook 
& Haehnelt 2006) do not rule out bubbles in the parameter range 
amenable for detection with the GMRT and MWA. 

In conclusion, we find that the optimal redshift for bubble de- 
tection is determined by a combination of instrument dependent 
factors and the evolution of the neutral fraction xhi. We propose 
that the redshift 8.1 ± 1.1 and 9.8 ± 1 are optimum for detect- 
ing ionized bubbles with the GMRT and MWA respectively. The 
prospects of a detection are most favourable for late reionization 
with xm ^ 0.5 at 2: 8 where for both instruments a 3a detection 
is possible for Rb > 30 Mpc with 1000 hrs of observation. Future 
observations will either lead to the detection of ionized bubbles, 
or lead to constraints on the product xni Rb for the observational 
voltune in the event of non-detection. 
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